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Abstract 

We study the effect of spatially dependent mass functions over the solution of the Klein-Gordon 
equation in the (3 + l)-dimensions for spinless bosonic particles where the mixed scalar- vector 
Coulomb-like field potentials and masses are directly proportional and inversely proportional to 
the distance from force center. The exact bound state energy eigenvalues and the corresponding 
wave functions of the Klein-Gordon equation for mixed scalar-vector and pure scalar Coulomb-like 
field potentials are obtained by means of the Nikiforov-Uvarov (NU) method. The energy spectrum 
is discussed for different scalar-vector potential mixing cases and also for constant mass case. 
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I. INTRODUCTION 



The relativistic wave equations can be considered as first approximation to the field theory 
when the corrections are taken in the presence of strong potential fields [1-3] . This explains 
the increased interest in the Klein-Gordon (KG) and Dirac wave equations to find exact and 
analytic solutions for energy spectrum and wavefunctions. As a consequence of the physical 
importance of exact solutions of relativistic KG wave equation in quantum mechanics, under 
the influence of strong potentials, an increasing interest in this equation has appeared in the 
study of quark- antiquark mass spectroscopy [4-6], atomic, nuclear, and plasma physics [7-9]. 
So, the idea is to use the bosonic particle KG equation as a mathematical tool to reach the 
goal of obtaining approximate solutions for nonrelativistic fermionic particles eigenstates. 
The problems that can be solved exactly in relativistic quantum mechanics are very limited 
perhaps because of the mathematical difficulties. 

Recently, the bound and scattering solutions of the s- and /-waves KG and Dirac equations 
for any interaction system have raised a great interest [10-12]. The bound-states of the Dirac 
and KG equations with the Coulomb-like scalar plus vector potentials have been studied in 
arbitrary dimension [13-17]. Furthermore, the exact results for the scattering states of 
the KG equation with Coulomb-like scalar plus vector potentials have been investigated 
in an arbitrary dimension [18]. This equation has been exactly solved for a larger class of 
linear, exponential and linear plus Coulomb potentials to determine the bound state energy 
spectrum using two semiclassical methods [19] . Many authors have considered a more general 
transformation between the unequal vector and scalar potentials given by 

V{r) = Vo + (5S{r), (1) 

where Vq and /3 being arbitrary constants of certain proportions have to be chosen after 
solving the problem under consideration [19-21]. It is interesting to note that, this restriction 
includes the case where V{r) — 0, when both constants vanish, the situation where the 
potentials are equal in magnitude and sign V{r) = S{r) or equal in magnitude but opposite 
in sign V{r) = —S{r) {i.e., Vq = 0; f3 = ±1), and also the case where the potentials are 
proportional when Vq = [21]. For the case when S{r) > V{r), there exist bound state 
solutions. However, it should be noticed that the case when the scalar potential is equal to 
the vector potential must be considered separately [3]. Under the condition of S{r) —V{r), 
the KG turns into a Schrodinger-like equation and thus the bound state solutions are very 
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easily to obtain with the help of the well-known methods developed in the non-relativistic 
quantum mechanics. 

On the other hand, the problem of the spatially-dependent effective mass is presenting a 
growing interest along the last few years [22-25] . Many authors have used different methods 
to study the partially exactly solvable and exactly solvable Schrodinger, KG and Dirac equa- 
tions in the presence of variable mass having a suitable mass distribution functions in ID, 3D 
and/or any D-dimensional cases for different potentials, such as the hnear potential [20], the 
exponential- type potentials [21], the Coulomb potential [26], the Lorentz scalar interactions 
[27], the hyperbohc-type potentials [28], the Morse potential [29], the Poschl-Teller poten- 
tial [30], the inversely linear scalar potential [31], the Coulomb and harmonic potentials [32], 
the modified Kratzcr-typc, rotationally corrected Morse potentials [33], Mic-type and pseu- 
doharmonic potentials [34]. Recently, the point canonical transformation (PCT) has also 
been employed to solve the £)-dimensional spatially dependent mass Schrodinger equation 
for some molecular potentials to get the exact bound state solutions including the energy 
spectrum and corresponding wave functions [32-34]. It is quite natural to look for relativistic 
treatment of this type of systems, mostly because the ordering ambiguity which is present 
in the nonrelativistic case [35], is expected to be avoided under relativistic ambiance [36,37]. 

Very recently, the NU method has been used to solve any i-states KG equation approxi- 
mately for the Hulthen potential with a suitable choice of spatially-dependent mass function 
distribution of an exponential- type [38]. Also, a new approximation scheme [39] has been 
proposed for the centrifugal term to obtain a quasi-exact analytic bound-state solution of the 
radial KG equation with spatially-dependent effective mass for scalar and vector Hulthen 
potentials in any arbitrary dimension and orbital angular momentum quantum number / 
within the framework of the NU method [40]. 

The problem of a particle subject to an inversely linear potential in one spatial dimension 
(~ known as the one-dimensional hydrogen atom, by considering a convenient mixing 

of vector and scalar Lorentz structures has received considerable attention in the literature 
[41,42]. The same problem for a spinless particle subject to a general mixing of vector 
and scalar inversely linear potentials in the (l+l)-dimensional world was analyzed [43]. 
Exact bounded solutions were found in closed form by imposing boundary conditions on the 
eigenfunctions which ensure that the effective Hamiltonian is Hermitian for all the points 
of the space. Jia and Souza Dutra [44] considered position-dependent effective mass Dirac 
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equations with PT and non- PT symmetric potentials. Souza Dutra and Jia [31], investigated 
the exact solution of the (l+l)-dimensional KG equation with spatially dependent mass for 
the inversely hnear potential. Recently, in ref. [45] the bound state solutions of the (1+1)- 
dimensional KG equation with mass inversely proportional to the distance from the force 
center for the inversely linear potential were obtained by using the NU method [46-52]. Two 
particular cases are studied, the case when vector potential is equal to the scalar potential 
in magnitude V{r) — S{r) and when vector potential is equal to the scalar potential in 
magnitude but not in sign V{r) — —S{r), (i.e., Vq — Q]I3 — ±1). 

In the present work, we feel tempted to extend the work of ref. [45] to study the bound 
state solutions of the (3 + l)-dimensional KG equation with position-dependent bosonic 
mass function m(r) = mo (1 + \obr~^) where r 7^ [45] for the attractive scalar potential 
S{r) — —Ksr~^ with Ks — hcQs being the coupling constant, taking into consideration the 
general mixings of scalar and vector Lorentz structure potential given in eq. (1). Firstly, 
this choice of mass function together with potential mixings is mostly suitable for modeling 
some physical systems hke the Kratzer-type potentials [48] . Secondly, the motivation for this 
choice is due to the nature of the dominating Coulombic field between the two interacting 
nuclei at short distances. Thirdly, this choice enables one to solve the KG equation easily 
and elegantly. The general mixings of potentials include: (i) V{r) — — K^r~^ and S{r) — 0, 
which represents a 7r~ meson in a Coulomb field, (ii) V{r) — and S{r) — —Ksr~^, which 
has no experimental evidence, (iii) V{r) — S{r) — —Kr~^, {ks — — k, Vq — 0, (3 — 1) 
which represents not only a KG particle in an equally mixed Lorentz scalar and Lorentz 
vector potentials but also a Dirac particle in the same potential mixture, where I = j + 1/2 
and the radial KG wave function represents the radial large-component of the Dirac spinor 
[53]. (iv) V{r) = —K,yr~^ + Ar and S{r) = representing a tt~ meson in a Coulomb 
field perturbed by a linear Lorentz vector interaction Ar. Also, we consider the effect of a 
spatially dependent mass of the linear form m(r) = mor/L [45] on the solution of the (3-1-1)- 
dimensional KG equation for the Lorentz vector and scalar potentials of the form V{r) — 
and S{r) = sr^^, respectively. It is worth mentioning that this choice of mass function 
together with potential mixings is mostly suitable for modeling some physical systems like 
the Pseudoharmonic potential. 

The paper is organized as follows. In sect. 2, we outline the NU method. Section 3 is 
devoted for the bound state analytic solutions of the (3 -|- 1) -dimensional KG equation with 
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spatially dependent mass functions for two quantum systems obtained by means of the NU 
method. Finally, the relevant results are discussed in sect. 4. 



II. NU METHOD 

The NU method is briefly outlined here and the details can be found in ref. [46]. This 
method is proposed to solve the second-order differential equation of the hypergeometric 
type: 

where a{z) and cr^z) are polynomials, at most, of second-degree, and r(s) is a flrst-degree 
polynomial. In order to find a particular solution for eq. (2), let us decompose the wave- 
function ipni^) as follows: 

'^ni^) ^<l>n{^)yn{z), (3) 

and use 

[a{z)p{z)]' = r{z)p{z), (4) 

to reduce eq. (2) to the form 

a{z)y:{z) + T{z)y'^{z) + Xy^iz) = 0, (5) 

with 

r{z) = t{z) + 27r{z), t'{z) < 0, (6) 

where the prime denotes the differentiation with respect to z. One is looking for a family of 
solutions corresponding to 

X^X„^-nT'{z)-^n{n-l)a"{z), n = 0, 1, 2, • • • , (7) 

The yn{z) can be expressed in terms of the Rodrigues relation: 

where Bn is the normalization constant and the weight function p{z) is the solution of the 
differential equation (4). The other part of the wavefunction (3) must satisfy the following 
logarithmic equation 

(9) 



a{z)' 



By defining 

k^X-7r'{z). (10) 

one obtains the polynomial 

Az) = \ W{z) - t{z)] ± ^\ [a'{z)-T{z)f-a{z)+ka{zl (11) 

where 7r(z) is a parameter at most of order 1. The expression under the square root sign in 
the above equation can be arranged as a polynomial of second order where its discriminant 
is zero. In this regard, an equation for k is being obtained. After solving such an equation, 
the k values are determined through the NU method. 

III. EXACT BOUND-STATE SOLUTIONS 

In the relativistic quantum mechanics, for a spinless particle, we write the full stationary 
KG equation for a spatially dependent bosonic mass in real (3 + l)-dimensions as [49,50] 

^'^(r) + ^ {[Em - V{T)f - [m{Ty + -S(r)]'} ^(r) = 0, = ^ ^ (12) 

j=l 3 

where m(r) is a bosonic mass, E^i is the energy of the particle, V{r) is a Lorentz vector 
(coupled as the 0-component of the four- vector potential) and S{r) is a Lorentz scalar (added 
to the mass term) potentials. Let us decompose the radial wavefunction ■0(r) as follows: 

^(r) = ^y«(?), (13) 

where u{r) is the radial wave function and Ym{v) is the angular dependent spherical 
harmonics and this reduces eq. (12) into the following position-dependent effective mass 
Schrodinger-like equation: 

+ ^ [eI, + V{rf - 2EV{r) - m{rfc' - S{rf - 2m{TyS{r) - 

X u[r) = 0, 14(0) = 0. (14) 

Now, we will start to analyze some illustrative particular cases for Lorentz scalar-vector 
mixings and suitable mass distribution functions for modeling some important physical sys- 
tems. 
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A. Mixed vector-scalar Coulomb potentials 



We shall present a general solution of ref. [46] for general admixture of scalar and vector 
potential mixings. Let us solve eq. (14) for the general relationship between vector and 
scalar potentials given in eq. (1), then we have 



(Pu{r) 1 



+ 



(1 - P') S{rf - 2 {{Er,i + Vo) /3 + m(r)) S{r) + 
1 



/(/ + l) ^V 



u{r) 



^2^2 [i^ni + V^Y - m^(r)c"] u{t). (15) 
Furthermore, we take the scalar potential in the form of an attractive Coulomb-like field 

%) = -^, = r^O (16) 
r 

where qs is being a scalar dimcnsionless real parameter coupling constant and he is be- 
ing a constant with J.fm dimension. At this stage, it is worthwhile to mention that the 
above choice of the mass function of Coulombic form together with the presently taken ad- 
mixture of scalar and vector fields are mostly suitable for modeling the well-known pseudo- 
Coulomb (Kratzer-type) potentials [48]. Following refs. [38,40,45], we may take the spatially- 
dependent mass function being described by 

m(r) = mo fl + ^ V r ^ 0, (17) 



r 

where mo and Aq = are the rest mass of the bosonic particle and the Compton-like 
wavelength in fm units, respectively; 6 is a dimcnsionless real constant. The interaction 
field has much impact on the choice of the mass function which, in the present case, is 
inveresely proportional to the distance between the two nuclei at short distances m(r) ~ ^ 
and constant at long distances m(r — >■ oo) ~ mo. Let us introduce the variable change z = r 
e (0, oo) and define 



^oc^ -Ki , 2 ^ X 2 (6 - g) moc2 - 2qj3E, 



Sni = 7:; {moc > E^i), 7i 



nl 



Q - Q 

Eni = Eni + Vo, 72 = M&-2g) + g' (1-/3') +/(/ + !), Q = he, (18) 

with the following constraint Vq < —Eni + moc^ must be fulfilled for bound state solutions. 
Further, substituting eqs. (16)-(18) into eq. (15), we obtain 

d^u{z) /4i^^ + 7i^ + 72~ 



( ^■"+J + "- )^M = 0. u(0)=0. (19) 



where u{z) — u{r). In the present work, we deal with bound state solutions, the quantum 
condition is obtained from the finiteness of the solution at infinity, i.e., the wave function 
u{r) must satisfy boundary conditions, u{r) =0 when r ^ oo and at the origin point, r = 0. 
In order to solve eq. (19) by means of the NU method, we should compare it with eq. (2). 
The following values for parameters are found: 

t{z) = 0, a{z) = z, a{z) = - [el.z' + + 72) • (20) 

Inserting these values of parameters into eq. (11), we obtain 

7r(z) l^Ki^^ + 4(7i + k)z + + 1. (21) 

The discriminant of the square root must be set equal to zero, i.e., A = '^s'l^iZ^ + 4(7^^ + 
k)z + 472 + 1 = 0. Consequently, the following two constants ki and k2 are found to be 

A;i,2 = -7i±£nzVl + 472, (22) 



with the following requirements on the parameters (3 < 2^^J {I + + [q — hY and Vq < 
—Eni + moc^ must be fulfilled for real solutions. In this regard, we can find the possible 
functions for n^z) as 

H ± [^niz + ix/1 + 472] for ki = -7i + £nl^/T+^^, 
t:{z) — < {Z6) 

y\± [sniz - 1^1 + 472] for k2 = -7i - SnlV^ + 472- 

According to the NU method, one of the four values of the polynomial 7r{z) is just proper 
to obtain the energy states because t{z) has a negative derivative for this value of 7r{z). 
Therefore, the selected forms of 7r{z) and k take the following particular values 

tt{z) = -Eniz + ^ (1 + a/I + 472) , A; = -7i - + (24) 

to obtain 

r{z) = -2eniz + 1 + ^/T+W2, r'(z) = -2eni < 0, (25) 

where t'{z) = In addition, after using cqs. (24) and (25) together with the assignments 
given in eq. (20), the following expressions for A are obtained 

Xn^X^2neni, n = 0,l,2,-- - , (26) 

A = -7i - eni (1 + a/T+I^) . (27) 



Letting A„ = A, we can solve the above equations for the energy states E^^ as 



q(b-q)(3± Br^i^Bl, - q^{l -(3^)-h{h- 2q) 



where 



Bni = n 



+ l + + +b{b-2q) + q^{l-p'), n, / = 0, 1, 2, • • • 



(28) 



(29) 



For spatially-dependent mass case, i.e., b ^ 0, the bound state solutions of the system are 
determined by the parameters q and b. It is not difficult to conclude that all bound-states 
appear in pairs, two energy solutions are valid for the particle — E^i and the second one 
corresponds to the anti-particle energy E"' — E'^ in the Coulomb-hke field. 

Let us now find the corresponding eigenfunctions for this system. Using eqs. (4) and (9), 
we find 

p{z) = ^\AT^e-'^"'^ (30) 
0(^) = ^(i+Vi^)/2e--"'^ (31) 
Hence, substituting eq. (30) into eq. (8), we find 



■„(z) = L>„z-v^e2^-^-^ L(n+v^)e-2^-^l ^ +i(2£„,z). 



(32) 



where L'^{x) is the generalized Laguerre polynomials. By using u{z) — (f){z)yn{z), we get 
the wavefunctions as 



u{r) = 7Vr(^+V'^)/2e--"<'-Lf +i(2£„^r). 



(33) 



where 



^{l+^^\b{b-2q) + q^ 



(34) 



Using the normalization condition u{r)'^dr = 1 and the orthogonality relation of the gen- 

r 1 2 

eralized Laguerre polynomials x°''^^e~^ lI^\x) dx — {2n -|- o; -|- 1) '"^""^""''"^^ , the nor- 
mahzing factor N can be found as [54-57] 



N 



n\ {2eni) 



2L+3 



2(n + L + l)r(n + 2L + 2) ' 
where Eni and L are given in eqs. (18) and (34), respectively. 



(35) 
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(1) If we consider the case when scalar potential is equal the vector potential in magnitude 
and sign, i.e., Vq — and P — 1, then we have 



rna 



q{b-q) + Bni^/Bl,-b{b-2q) 



q^ + B, 



nl 



q{b-q)- Bni^Bl-b{b-2q) 



g2 + 5,2 



moc^ 



moc^, 



nl 



where 



1 + +b{b-2q), n,/ = 0,l,2, 



(36a) 



(36b) 



(37) 



Obviously, the bound state solutions of the particle and anti-particle are available. When 
the mass is taken to be constant, i.e., 6 = 0, we have 



{n + l + lf-q^ 



mQ(?, E^i — —mQ(?, n,l — 0, 1, 2, 



(38) 



(n + l + iy + q^ 

where n and / signify the usual radial and orbital quantum numbers. The particle has bound 
state solution whereas anti-particle has continuum solution for all states. 

(2) If we consider the case when scalar potential is equal to the vector potential in 
magnitude but not in sign, i.e., Vq = and l3 — —1, then we have 



-q {b-q) + B^i^Bl,-b{b-2q) 



-q (b-q)- B^i^Bl,-b{b-2q) 



q^ + B"" 



nl 



For the constant- mass case, i.e., b — 0, we have 



n + l + iy-q^ 2 
moc . 



(39a) 



(39b) 



(40) 



{n + l + iy + q^ 

Obviously, the particle has continuum solution for all states whereas bound state solution 
for anti-particle. In addition, when the potential coupling constant is taken as q — b/2, the 
spectra of the varying mass KG particle in potential fields qs — qv {qs — —qv) are similar 
to the spectra of constant mass KG particle in the potential fields = —q^ (g^ = q^), 
respectively. 



10 



B. Pure Scalar Coulomb- like potential 



In their paper [44], Souza Dutra and Jia used the pure scalar potential that is inversely 
proportional to the absolute value of the coordinate. Here, we use a pure scalar repulsive 
Coulomb-hke field potential 

S{r) = V{r) = 0, (41) 

with s being a coupling parameter with J.fm dimension and also assume the spatially- 
dependent mass function having a linear form 

m(r) = Ar, (42) 

with A = where mo is being the rest mass and L is being a constant with space dimension. 
At this stage, it is worthwhile to mention that the above choice of the mass function together 
with the presently taken pure scalar potential case are mostly suitable for modeling the well- 
known pseudoharmonic potentials [58,59]. Inserting eqs. (41) and (42) into eq. (15), then 
we have 



d'^u(r) 1 



+ 



<r) -^AEli- ^ ) u{r). (43) 



12 ^2 

Thus, the present problem has been reduced to three-dimensional Schrodinger equation for 
pseudoharmonic oscillator problem which was solved before in refs. [58,59]. Introducing the 
variable change z — r'^ & (0, oo) and defining 



we obtain 



d'^uiz) 1 duiz) 1/22 2 \/\^ 

^ + T^-ir + (-"'^ - - "<^> = »■ 

where u{z) — u{r). Comparing eq. (45) with eq. (2), we find values for the parameters as 

t{z) — 1, a{z) — 2z^ 5(2;) = —a\z^ — e'^iZ — 0:2, (46) 

and by inserting these values of parameters into eq. (11), we further obtain 

1 , 1 



n{z) = -± -^4alz^ + 4(4 + k)z + 4^2 + 1, (47) 
and the constant k as 

ki,2 = -ell ± q;iV4q;2 + 1. (48) 
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When the individual values of k given in eq. (48) are being substituted into Eq. (47), the 
four possible forms of Tr{z) are written as follows 

I ± I [2aiz + V4q;2 + l] for ki = -eli + q;iV4q;2 + 1, 
■7r{z) = ^ (49) 

i ± I [2aiz - V4q;2 + l] for k2 = -eli - ctiViaT+T. 

According to the NU method, the selected forms of t:{z) and k are taking the following 
particular values 

n{z) = -aiz + ^ [l + \/4q;2 + Ij , k = -eli - q;iV4q!2 + 1, (50) 

to obtain 



t(^) = -2ai;2 + 2 + V4a2 + 1, t'(;s) = -2ai < 0, (51) 

which is the essential condition in the method. Also, the following expressions for A are 
obtained 

A„ = A = 2nai, n = 0, 1,2,--- , (52) 
A = -ell - ai (1 + V4q;2 + l) . (53) 
Letting A„ = A, we can solve the above equations for the energy eigenvalues Eni as 



E^i = ai (1 - 2A^i) , (54) 



where 



and thus we find 



E^i = -4, A^i = -{n + lJ{2l + lf + ^], (55) 



2 V ' ' K^c 



+ 2n+l + ,(2,+ lf + ifil. (56) 



moc^ L L \ V h?c^ ' ' 

which is found to be consistent with eq. (20) of ref. [51] obtained by SUSY method when 
I is set equal to zero. We note that the energy levels for particles and antiparticles are 
symmetric about Eni = [49]. 

Essentially, we should report that eq. (43) corresponds to the Schrodinger equation of 
anharmonic oscillator potential V{r) = afr^, with energy levels [58,59] 

E,; = ai(2n + 2A + 3), (57) 
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with Eni is given in (55), ai in (44) and A is defined by 




(58) 



Let us now find the corresponding eigenfunctions for this system. After using eqs. (4) and 
(9), we find 

6(z) = ^(l+V45^+T)/4g-ai./2_ 



Substituting eq. (59) into eq. (8), we obtain 



yn{z) = D^^-V4a2+l/2gai. 



(n+V4a2+l/2) „-aiz 



L(fA+i)/2(ci^). 



(60) 



(61) 



By using u{z) = (f){z)yn{z) , we get the wavefunctions as 

u{r) = Are-"i'-V2^(A+l)/2^(2A+l)/2(^^^2)^ 



(62) 



where A is defined in eq. (58). It is worth mentioning that the above wave function is 
consistent with eq. (20) of ref. [58] in the solution of the Schrodinger equation for the 
pseudoharmonic oscillator potential. Essentially, such a solution has been discussed before 
by many authors [59]. Making use of the normalization condition u{r)^dr = 1 and the 
orthogonality relation of the generalized Laguerre polynomials x 



dx 



a — a' + n 



!(«'), the normalization constant N can be found as [54-57] 



n 



N 



(63) 



n — 1 



n 



rU\/(2^ + i)+(S) +1 



IV. CONCLUSIONS 

We have extended the (l+l)-dimensional KG solution in ref. [45] to the i-waves KG for 
scalar-vector mixing Coulomb-like fields with suitable choices of spatially dependent mass 
functions. Thus, for this kind of studied problems, we may conclude that the relativistic 
wave equation can be solved exactly. For suitable choices of potential forms as the general 
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mixing of scalar- vector and pure scalar Coulomb-like field potential, the relativistic bound 
state energy spectrum and wave functions have been obtained, respectively. The resulting 
solutions of the wave functions are being expressed in terms of the generalized Laguerre 
polynomials. We have considered different mass functions of inversely proportional and 
directly proportional to the coordinate distance. Obviously, when the coupling potential 
parameters are adjusted to some specific values, particularily when q — h/2, the spectra of 
the mass varying KG particle for the case qs — Qv {Qs — —Qv) become similar to the spectra of 
the constant mass KG particle for the case Qg — —Qv {Qs — Qv), respectively. It is found that 
the KG equation with a suitable mass function for a pure scalar potential is being reduced 
into the constant mass Schrodinger equation for the anharmonic oscillator potential. In the 
limit of constant mass (6 = 0), the solution for the energy eigenvalues and wave functions 
are reduced to those ones in literature. Also, when I — 0, the problem reduces to s- waves 
solution. 
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